Probing spatial spin correlations of ultracold gases by quantum noise spectroscopy 
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Spin noise spectroscopy with a single laser beam is demonstrated theoretically to provide a direct 
probe of the spatial correlations of cold fermionic gases. We show how the generic many-body 
phenomena of anti-bunching, pairing, antiferromagnetic, and algebraic spin liquid correlations can 
be revealed by measuring the spin noise as a function of laser width, temperature, and frequency. 
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Ultracold atoms offer the possibility to prepare, manip- 
ulate and probe various paradigm phases of strongly cor- 
related systems. Considerable efforts are devoted to de- 
velop sensitive detection schemes to study these phases. 
Whereas most experiments in this field are based on mea- 
suring mean values of various observables, further insight 
can be obtained from the correlations in the noise of the 
atomic distribution [J, 0, Q ■ recent experiments a new 
technique using phase contrast imaging was used to probe 
the spin of ultracold atoms [1, Q . In related experiments 
d, 0,11] similar techniques have been pushed to the point 
where they are sensitive to the quantum fluctuations of 
the atoms. In this Letter, we show that quantum spin 
noise spectroscopy along the lines of Refs. 0, H, 0, 
constitutes a sensitive probe of the correlations of the un- 
derlying quantum state. We focus on generic many-body 
phenomena such as antibunching, pairing, and spin liq- 
uids. Furthermore, we show that spin noise measurement 
is an ideal tool for probing antiferromagnetic ordering 
and phase transitions for atoms in optical lattices, which 
is currently one of the main challenges in cold atoms 
physics. Related theoretical studies of spin noise have 
recently been presented in Refs. [ol. [lol. . 



Quantum noise limited probing of the spin state may 
be obtained either by polarization rotation 01 or phase 
contrast imaging 0, . In the first approach the spin 
imprints a phase shift on a laser beam and this phase shift 
is subsequently measured by interfering the beam with 
another laser beam (i.e. homodyne detection). In the 
polarization rotation measurement the two laser beams 
are replaced by two different polarization modes, which 
has the advantage that the setup is less sensitive to fluc- 
tuations in optical path length and beam proflle. Ideally 
one would probe the system by imaging with a camera, 
but here we explore a slightly simpler situation where a 
laser beam is passed through the sample and the flnal re- 
sult is measured by photo detectors without any spatial 
resolution. In the limit of strong beams (many photons) 
experiencing a small phase change, the observable, i.e. 
the measured light quadrature in the homodyne detec- 



tion, may be expressed as 0, [n, M, m 
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Here, Xin/out is a canonical position operator describ- 
ing the light normalized such that the input corresponds 
to vacuum noise (X?^) = 1/2, and k is a coupling 
constant. The effective measured atomic operator is 
= J (fr (t){r)szir)/VA, where A = / d^r(?!)2(r)?i(r)/4 
is a normalization constant, (j){r) is the spatial intensity 
proflle of the laser beam, and Sz(r) = (?/'j(r)-0| (r) — 
'0|(r)'0|(r))/2 ~ {n^{r) — 7i|(r))/2 gives the local popu- 
lation imbalance (magnetization) with if^ai^) being the 
atomic fleld operator. We consider a two-component 
atomic gas (cr =t,i) with total local density n(r) = 
(ri|(r)-|-7i|(r)) and assume Gaussian laser profiles 0(r) oc 
g-(a: +y )/d ^ gy measuring the observable Xont it is pos- 
sible to obtain spatially resolved information about the 
magnetization (Xout) = k,{Mz) / V^. In many cases, how- 
ever, interesting states may not have any net magnetiza- 
tion (Mz) = 0. In this Letter, we will only consider such 
situations and show that a measurement of the quantum 
(1 + Rk^)/2, where R = (M^ 
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provides insight into the state of the system. Since R is 
quadratic in the atomic density operators it gives a direct 
measure of the atomic correlations in the system. The 
normalization in Eq. ([2]) is chosen such that the quantum 
noise of an uncorrelated state, where each atom has an 
equal probability of being in each of the two internal 
states, is i? = 1 (standard quantum limit). 

We first consider the normal phase, where the spin fluc- 
tuations have a length scale of fc^^ . It follows that R van- 
ishes if the effective volume Vb — {Jy cPrcj))'^ / Jyd^rcf)'^ 
is large, Vb ^ kp^. Fermi statistics thus suppresses the 
noise below the standard quantum limit R = 1. For a 
finite laser beam there will, however, be a noise contri- 
bution from the boundary (A/_^) ^ d, which translate 
into R ^ l/kpd. 
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A key property of pairing for fermions is that the two 

particle density matrix (^/'|(ri)V'|(r2)'0i(r'2)^T('''i)) ^'"^^ 
a macroscopic eigenvalue PcN with N the number of par- 
ticles and Pc the condensate fraction. Spin noise spec- 
troscopy probes the two particle density matrix directly, 
and in the large d limit the noise is dominated by the 
largest eigenvalue PcN. The noise depends on the shape 
of the applied laser beam as seen from the following ar- 
gument: assuming a top hat laser profile with a radius 
d and sharp edges compared to the radius ^ (coherence 
length) of the pair wavefunction x(i")i the noise is pro- 
portional to the number of pairs within ^ of the edge 
such that only one particle is inside the beam. This 
gives a scaling R oc 1/d as in the normal case. With a 
smooth laser profile with radius d and fall-off distance 
D > ^, the noise is due to pairs in the edge region. 
These pairs couple to the gradient l/f) and the noise 
from the difference in signal from f and j particles is 
~ / d^rx^{r)r'^ /D'^ - ^V-D^- This should be multiphed 
by the number of pairs in the edge region ~ LzdDpcN/V, 
where and V denote the length and volume of the sys- 
tem. Since A ~ L^d'^N/V, we get R ~ pd'^/Dd. With a 
Gaussian beam D ^ d, and the scaling R pc^"^ /d"^ thus 
provides a measurement of ^ and pc- 

We now use the BCS wavefunction to derive this scal- 
ing rigorously in the BCS and BEC limits. Consider a 
homogeneous gas with constant density ncr(r) = N„/V . 
Wicks theorem yields {sz{yi)sz{y2)) = n(5(r) — 29^{r) — 
2F^{r) with r = n - ra, 6'(r) (vSj (n (ra)), and 
F{r) = (V)|(ri)V;^(r2)). We then find 

R=l-jJ dVid3^20(ri)0(r2)[02(r) + F^{r)]. (3) 

In the BEC regime kpa 0+, the chemical potential is 
^ — ft^/2ma^. This gives UkVk A/2(|/i| -I- fc^/2m) 
and w| A^/4(|/i| -I- /2mY for the coherence factors 
defined as = (1 + ^/E)/2, v'^ = 1 - u"^ with E = 
(^2 + A2)i/2 and ^ k^/2m - n We ob tain e{r)/n^ = 
exp(— r/a) and F(r)lna ~ ^Ji-n/kpa eyip{— r /a)/ kpr 
which is proportional to the asymptotic bound state 
wavefunction for a potential with scattering length a. 
Likewise, in the BCS limit kpa — > 0_, 6{r)/n„ = 
3[sinA:F'' — kpr cos kpr^ irr /2^ exp(— r/^)]/(fci?r)'^ and 
F{r)/n^ = 3sinfci?r-y/7r72^exp(— r/^)/fc|,r for r oo 
where ^ = kp/niA (h = 1) and A is the gap. Using these 
limiting forms in ([3]), we obtain for d ^ oo 
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For s-wave interactions, the pair wavefunction has a 
short-range divergence (bunching) given by F{r) = 
mA/47rr [15| resulting in a linear decrease of the noise 
for fci?(i — > in both the BCS and BEC limits. Using 
Pc ~ l/kp£^, the BCS result agrees with the estimate 
given in the previous section. 



1 



oc 



o. 

















» 'i \ " -. \ \ 0.5 
















' \ \ .\ \ O 


SO 1 OO 






Normal 

k^a=-1 
















k^a=0.5 













o 1 2 ^kd'* ^ ^ 

FIG. 1: (Color online) Noise R{d) for various coupling 
strengths. The thin lines show the d ^ oo limit @. The 
inset shows the large d behavior d/L ^ Oil) with P = 1 for 
the normal phase and the superfluid phase in the BEC limit. 

In Fig. [1] we plot R{d) for a homogeneous system 
of transverse radius L ~ lOO/s^^^ at T = 0. Results 
for the normal phase and the superfluid phase with 
{kpo)^^ = —1 (BCS regime), {kpa)~^ = (unitary 
limit), and {kpo)^^ = 2 (BEC regime) are shown. The 
noise is calculated numerically from ^ using the BCS 
wavefunction. The noise is below the quantum limit 
i? < 1 and i? ^ for d ^ oo in agreement with the 
analysis above. Pairing suppresses the noise compared 
to the normal state due to positive correlations between 
opposite spin states. The suppression increases with the 
pairing moving toward the BEC side. 

For very large d the laser probes a significant frac- 
tion of the system and it is important to include pos- 
sible spin fluctuations due to the experimental prepa- 
ration of the system. Typically, such fluctuations will 
at least be limited by the standard quantum limit, i.e. 
R — P with (3 > 1 when probing the entire system. 
When probing a sub-system this gives an extra contri- 
bution ^ (3Vb/V which is important for large d for the 
normal phase and the superfluid phase on the BEC side 
{kpa)~^ = 2 (see inset in Fig. [1]). However, this term is 
absent for {kpa)~^ < 0.5, since superfluidity quenches 
the spin noise in this regime 

U B El- ■ 

Observing 

R 1 for a large portion of the sample would represent 
an extreme experimental demonstration of this quench- 
ing. 

The observed enhancement of the nuclear spin relax- 
ation just below the transition temperature Tc (Hebel- 
Slichter effect) constitutes one of the hallmark experi- 
mental tests of BCS theory. We now demonstrate the ex- 
istence of a spin noise spectroscopy analogy to the Hebel- 
Slichter effect. Similar effects has been demonstrated to 
occur in inelastic light scattering and Bragg scattering 
experiments [31 . The probing technique discussed in 
this Letter is in principle non-destructive. By recoding 
the signal for a long duration of time one can thus obtain 
all frequency components of the noise R{d, to) by Fourier 
analysis [lO|, i.e. Fourier transforming the measured 
Xout(t) provides a measurement of Mz{uj). Such probing 
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FIG. 2: (Color online) R{d,uj) in units of R{d,uj) at Tc for 
various laser widths d and kpa = —1. 



will have similar signal-to-noise ratio 



but 



since spontaneous emission may lead to significant heat- 
ing 77 may have to be kept very low to avoid that the 
system heats up during the measurement. Using ([3]) we 
obtain for a homogeneous system 
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where Iq is the modified Bessel function of the first 
kind, kx — (k^yky) is the transverse momentum, and 
/ = [exp{(3E) + 1]^^. The primed quantities refer to the 
momentum k' with energy E' = E + lo. There is mo- 
mentum conservation along the z-direction with k'^ = k^ 
whereas k'^ ^ kj^ due to the transverse Gaussian profile. 
Eq. ([5]) gives the noise contribution from quasiparticle 
scattering from momentum k to k'. There are additional 
terms describing pair breaking and quasiparticle absorp- 
tion which do not affect the Hebel-Slichter effect. 

In Fig. [21 we plot R{d,uj) as a function of T/Tc cal- 
culated numerically from ([5]) using the self-consistently 
determined gap A(T). We have chosen kpa = —1 giv- 
ing Tc/Tf ~ 0.13 and w/fcsTc = 0.08 since the Hebel- 
Slichtcr effect only occurs for lu < ksTc. For narrow laser 
widths, a Hebcl-Slichter peak is prominent below Tc- The 
peak arises from an increased density of states at the gap 
edge; it decreases with increasing d and disappears for 
d » C(T = 0) ~ 9. This is because for large laser widths, 
the scattering becomes subject to momentum conserva- 
tion which restricts the available phase space. 

Quantum systems in periodic potentials constitute an- 
other class of intriguing systems which can be examined 
by cold atomic gases using optical lattices. Superfluid- 
ity in lattices, possibly of d-wave symmetry, can be de- 
tected by suppression of spin noise similar to the discus- 
sion above for homogeneous systems. The only difference 
is that the d-wave pair wavefunction does not diverge for 
short length scales and there is no linear decrease in R for 
small laser radii d. One could use a laser with elliptical 
transverse profile to detect the anisotropic suppression of 
spin noise due to the d-wave symmetry of the pairing. 



Presently, a main experimental goal in optical lattices 
is to observe the onset of antiferrornagnetic (AF) correla- 
tions with decreasing temperature [20] . As demonstrated 
below, spin noise spectroscopy can measure the magnetic 
susceptibility of the system and hence constitutes an im- 
portant experimental probe of the spin correlations. As 
an example, we study atoms described by the Hubbard 
model which in the strong repulsion limit at half fill- 
ing for kT ^ U reduces to the AF Heisenberg model, 
H = J j) • Sj, where {i, j) denotes nearest neighbor 
pairs and is the spin 1/2 operator for the atoms at 
site i. Assuming, without loss of generality, a staggered 
magnetization along the z-direction, we now show how 
to detect AF correlations by measuring = (M^) and 
R± = (A/J), where is defined analogous to M^. (A 
preferred direction for the broken symmetry can be in- 
duced by enforcing a slight anisotropy in the exchange 
coupling J.) Here we are mainly interested in the T de- 
pendence and focus on the situation where we probe the 
entire ensemble. Therefore, we assume a broad laser pro- 



1 in m such that R 



{x)Sz{x))/N 



file with (j: 

with S = Si and N is the number of spins. In the 
paramagnetic phase, R\\ = R± ~ AkTx where x is the 
magnetic susceptibility. A high temperature expansion 
yields for the 2D square and 3D cubic lattices 21 1 



AkTx = 



1 - 2.x 
l-3x 



2x2 _ 1 333a;3 



2D 
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where x ~ J/2kT. In 2D, the system remains param- 
agnetic for T > 0, and modified spin- wave theory yields 
X = (12J)-i[0.524-h0.475r/J-He'(T3)] for T/J < 1 
In the 3D case, the system undergoes a phase transition 
to an AF phase at the Neel temperature T/y. In the AF 
phase, R±. Using spin- wave theory for T < Tj^, 

we obtain R± = kT/{3J) and 



Rn 



1 
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N ^ 2sinh^(/3wk/2)' 



(7) 



Here = 3 J-\/l — 7^ is the spin- wave energy with 
7k = [co^kr^a -I- cos ky a + cosfcza)/3 for a cubic lat- 
tice with lattice constant a. The sum in ([7|) is over 
the reduced Brillouin zone. For kT <C J, (O yields 

= 4(A:T)3/(3s3), with s = V3J the spin-wave veloc- 
ity. In Fig.[3lwc plot these results for both the 2D and 3D 
systems. We see that the onset of AF correlations in the 
paramagnetic phase can be detected as a decrease in the 
noise from the uncorrclatcd result _R|| = R± = 1 as de- 
scribed by By comparing with the high temperature 
expansion, the spin noise may even serve as an accurate 
thermometer for the spin temperature. Furthermore, the 
AF phase for the 3D case can be detected by observing 

^ R±. An advantage of probing collective operators 
like Sz is that they are conserved, and therefore could 
be measured after time of flight. In this case, however. 
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FIG. 3: (Color online) and Ri_ for a 2D (a) and 3D (b) 
system. Solid lines are the high and low T results discussed 
in the text and the dashed lines in (b) (obtained by simple 
rescaling) indicate how they meet at Tjv — 0.946 J [2a |. 



special care has to be taken of the contribution from the 
boundary region. We do not expect the trapping poten- 
tial to change these results quaUtatively 



M 2 



The method presented here can also be used to probe 
the correlations of more exotic quantum phases such as 
resonating valence bond states and algebraic spin liquids 
[2^. These states are characterized by long-range spin 
correlations (s(r)s(O)) ~ (— l)''^+''y/r(^+''). Techniques 
exist for addressing, e.g., every second site in a lattice 
(27| . Flipping every second spin before the measurement 
(s^ ^ (-!)'■- will give (s^(r)s^(O)) ~ l/r(i+''). Per- 
forming noise spectroscopy on this state will give a contri- 
bution from the long-range correlations R ~ d^^~^\ By 
measuring the scaling of R with d one can thus directly 
determine the exponent 77 of the spin correlations. 

Finally we consider the experimental requirements for 
realizing our scheme. The experiments should be quan- 
tum noise limited with all classical noise sources sup- 
pressed. This has already been achieved in several ex- 
periments [1, 0, S, , and we expect it to be simpler to 
realize for the smaller systems considered here. In addi- 
tion, the atomic noise should be large compared to the 
light noise inherently present in the probe. The sponta- 
neous emission probability pr. atom caused by the prob- 
ing light is 77 ^ K^/a, where a = 3nL^A^7a;/727r is the 
optical depth of the ensemble 0,[i3,[il. Taking n - 10^^ 
cm~'^, Lz ^ 100 /im, a probing wavelength A ~ 671 nm 
corresponding to Li, and a branching ratio "fx/^ ~ 1/2 
gives a = 16 for a harmonically trapped Fermi gas. For 
atoms in optical lattices at half filling a w iV^ = 50 where 
Ns is the number of lattice sites in each direction. One 
can thus have a large signal-to-noise ratio k^R> 1 with 
very little noise added from spontaneous emission dur- 
ing the probing 77 = n^/a <C 1. Another concern is the 
spatial resolution. Experimentally, one may obtain a res- 
olution down to d ~ 5A Taking n ~ lO^^cm^"^ this 
corresponds to kpd ~ 10. Thus, it may require an adia- 
batic expansion of the gas to observe the small scale limit 



of Fig. [T] However, it is possible directly to observe the 
large d scaling, the Hebel-Slichter effect, and the onset of 
magnetic correlations. 

In summary, we have shown how to extract the corre- 
lations of quantum states of ultracold atoms using spin 
noise spectroscopy. This was demonstrated explicitly by 
calculating the spin noise for normal Fermi gasses, super- 
fluids, paramagnetic and AF phases and algebraic spin 
liquids. This method can be applied to other strongly 
correlated systems as well as extended to higher order 
moments [ll| . It may even be extended to full quantum 
state tomography of the two particle density matrix. 
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